Localization of strongly correlated electrons as Jahn- Teller polarons in manganites 
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A realistic modeling of manganites should include the Coulomb repulsion between Cg electrons, the 
Hund's rule coupling to t2g spins, and Jahn- Teller phonons. Solving such a model by dynamical mean 
field theory, we report large magnetoresistances and spectra in good agreement with experiments. 
The physics of the unusual, insulating-like paramagnetic phase is determined by correlated electrons 
which are — due to strong correlations — easily trapped as Jahn- Teller polarons. 

PACS numbers: 71.30.-|-h, 71.27.-|-a, 75.47.Gk 
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Since large magnetoresistances are of substantial tech- 
nological importance, the colossal magnetoresistance 
(CMR) Q in manganites (Ti-^D^jMnOa; T: trivalent 
ion, e.g., La; D: divalent ion, e.g., Ca) has attracted 
considerable attention, both experimentally and theoret- 
ically Like the mechanism for high-temperature su- 
perconductivity in the cuprates, the CMR still lacks a 
thorough theoretical understanding; and like supercon- 
ducting transition temperatures, the temperatures for 
the CMR are much lower than technologically desirable. 

The key to the understanding of CMR lies in the 
unusual properties of the paramagnetic (PM) phase 
which — in a wide range of dopings x — shows an unusual 
insulating-like behavior. At lower temperatures T and/or 
within magnetic fields, the double exchange mechanism 
stabilizes a ferromagnetic (FM) metallic phase. If 
the insulating-like PM phase was understood, a CMR 
could simply be described as the PM-insulator-to-FM- 
metal transition. Different proposals aim at explaining 
the CMR and the PM insulating-like phase, in partic- 
ular, the localization of charge carriers through lattice 
or orbital (5iJ polarons, an Anderson localization aris- 
ing from disordered i2g-spins |^ , a phase separation into 
nano-domains showing percolation effects 0, [1] , and an 
effective two band model However, hitherto no quan- 
titative, microscopic calculation satisfactorily explains 
the known physical properties of the PM insulating-like 
phase, which is not only characterized by an increase of 
the resistivity p with decreasing T, but also has unusual 
dynamic properties. The latter reflect in a spectral func- 
tion A((jj) with a very low spectral weight at the Fermi 
level Ep irrespectively of indicated by photoemis- 
sion and X-ray absorption experiments [lolllll| . Similarly 
the optical conductivity a{uj) shows a very low spectral 
weight up to an energy scale of ~ 1 eV . Also the FM 
is an atypical (bad) metal. 

The theoretical understanding is complicated since 
several ingredients seem to be necessary to describe man- 
ganites: the exchange interaction between more local- 
ized t2g spins and the itinerant eg electrons 0, the 
Jahn- Teller phonons and their coupling to the electrons 
0, as well as the electronic correlations due to the lo- 
cal Coulomb interaction 0, Q|. The necessity of all 



three of these interactions was also revealed by realistic 
LDA-(-DMFT (local density approximation -I- dynamical 
mean field theory |l5l |) calculations for the parent com- 
pound LaMnOa, showing the metal-insulator transition 
in LaMnOs under pressure 0| to be triggered by the 
(Jahn- Teller) crystal field splitting which is strongly en- 
hanced by the local Coulomb interaction [l7j |. 

In this Letter, we model doped LaMnOa by including a 
realistic tight- binding band structure for cubic LaMnOs, 
the Coulomb interaction between the Cg electrons, Jahn- 
Teller phonons, and the Hund's exchange coupling to lo- 
calized t2g spins. We solve this model by dynamical mean 
field theory (DMFT) employing the (numerically) 
exact quantum Monte Carlo (QMC) method as an 
impurity solver. We find that tendencies of the Coulomb 
repulsion and the Jahn- Teller distortion to localize elec- 
trons mutually support each other, trapping Cg electrons 
as local polarons in the PM phase. Hence the PM phase 
is insulating-like with very low optical spectral weight 
below leV and with a resistivity which is by a factor 
of ~ 8 larger than that of the FM phase (magnetoresis- 
tance). Jahn- Teller polarons without correlation effects 
have been studied before, e.g. in 0,113,01, but suffered 
from shortcomings such as big magnetoresistances only 
for undoped LaMnOa (which is always insulating in ex- 
periments) and discrepancies with experimental spectra. 

Model for manganites. Starting point of our investiga- 
tion is the following, realistic model for manganites 
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Here, cj^^. and c^^^ are creation and annihilation op- 
erators for electrons on site i within eg orbital m and 
spin cr; s,^ = hUaa' clua'^^^'Ci^^' dcuotes the e^-spin 



(r: Pauli matrices), Si the t2g-spin, and Qa {Pa) the 
coordinate (momentum) of the two quantum Jahn- Teller 
phonons. Let us briefly discuss the three lines of Hamilto- 
nian The first line describes the Kondo lattice model 
consisting of two terms: (i) the tight binding band struc- 
ture with directionally dependent nearest-neighbor hop- 
ping, tfi = -V3tf2 = -\/3tii=3tf2 = 3to/4, tl^^V3t\^ = 
VSt^i = 3t|2 = Sto/i, t^2 = to, and tf^ = tf^ = tfi = 
with bandwidth M^ = 6to = 3.6eV which well describes 
the LDA band structure for the cubic lattice (contribu- 
tions of longer-range hopping are very minor |22j|'l: and 
(ii) the coupling J' to the t2g spin which we assume to 
be classical with strength 2j7|Si| = 2.66 calculated from 
a ferromagnetic LDA calculation jT^I • The Kondo lattice 
model (first line) can be solved exactly in DMFT |2^ 
and gives rise to the double exchange mechanism for 
ferromagnetism. The second line describes the Coulomb 
interactions between the Cg electrons including the intra- 
([/) and inter-orbital repulsion [/', as well as the Hund's 
exchange (J). Even for large the Coulomb repul- 
sion U' — J between two spin-aligned Cg electrons on the 
same lattice site is important. It leads to the formation of 
Hubbard bands and quasiparticle peaks in the PM phase 

3. We take U' = 3.5 eV and J = 0.75 eV from Ref. 

nf (C/ = [/' -I- 2J follows by symmetry). Finally, the 
third hne describes the couphng g of the Cg electrons to 
the two local Jahn- Teller phonons; and the fourth line 
the quantum motion of these phonons. The phonon fre- 
quency — 0.07 eV is estimated from Raman spectra 
and the results of lattice dynamical calculations |2J| ; the 
breathing mode, which has a higher frequency and is be- 
lieved to be less important is neglected. Since the 
electron-phonon coupling is difficult to determine from 
band structure data, we consider different values, mostly 
(7 = 0.10eV'^/^ which corresponds to a dimensionless cou- 
pling constant A — g/ft^/t = 1.84, comparable — slightly 
larger — than previous estimates, see |8|. In the follow- 
ing, we set h = e = kB = ^', our unit of energy is eV. 

DMFT( QMC) implementation. For solving Hamilto- 
nian (HJ by DMFT(QMC), the imaginary time r = 0../3(ee 
1/T) is discretized into 1 — 1..L Trotter slices n of width 
Ar = 0.2, and the Coulomb interactions are decoupled 
through discrete Hubbard-Stratonovich transformations. 
We sample the auxiliary Hubbard-Stratonovich field and 
the continuous phonon Bose field Qa{Ti) on equal foot- 
ing, after replacing the phonon momentum Pa{Ti) by 
[Qain) — Qaiji-i)]/ /S.T as discussed in [i^. After each 
QMC sweep, a global update of the phonon field is con- 
sidered; we use ~10^— 10'' sweeps. 

Polaron formation in the PM phase. Let us start our 
discussion with the QMC distribution of the r-averaged 
phonon field Q2, corresponding to the amplitude of the 
lattice distortion [2^. Fig. ^ shows three peaks with 
large (positive and negative) Jahn- Teller distortion and 
small Jahn- Teller distortion, respectively. Upon lower- 
ing temperature, and at strong enough electron-phonon 
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FIG. 1: (Color online) QMC distribution of the lattice distor- 
tion [1/1/ X^j^i Q2{ti)\ for different electron-phonon couplings 
g, n = 0.8 electrons/site and /3 = 16; the dotted line is without 
Coulomb repulsion for comparion [g — Q.l). At strong enough 
g, three regions with large (positive and negative) and small 
lattice distortion emerge. The inset shows the electron spec- 
tral function for these regions with large (light/green) and 
small (dark/red) lattice distortions, respectively. We can in- 
terpret the results physically as the trapping of the n = 0.8 
electrons as Jahn- Teller polarons. 



coupling, these three peaks get sharper and clearly sepa- 
rated in the PM phase. Then, the QMC sampling tends 
to stay within one of the three regions for a very long 
computational time [2^ . 

Physically, the three sharp peaks mean that the local 
lattice site is either strongly Jahn- Teller distorted (with 
orbital 1 or 2 lower in energy as visualized in Fig.^l or not 
(central peak). This clear separation into three peaks fur- 
thermore suggests that the phonon dynamics is strongly 
reduced. In our view, the slow QMC dynamics for going 
from one of the two distorted to the undistorted config- 
urations indeed reflects a slow real time dynamics j28| . 
Turning to the electronic spectrum in the inset of Fig. ^ 
we see there is one electron with Jahn- Teller distortion 
(the peak below Ep = of the light/green line contains 
one electron) and almost none without distortion (the 
dark/red line is above the Fermi energy). This is ex- 
actly the picture of a local polaron, an electron trapped 
through the lattice distortion. This trapping explains 
why we obtain an insulating-like spectrum at g = Q.\ and 
even a true gap for 3 = 0.15. Let us note that the polaron 
formation is strongly supported by the local Coulomb re- 
pulsion. Without U' , the polaron formation fades away 
(dotted line in Fig. The effect of U' on the polaron 
formation is two-fold: (i) It strongly enhances the Jahn- 
Teller splitting from 2EjT=g^/^'^ to 2Ejt + U' - J, see 
the (light/green) spectrum in the presence of distortion 
in the inset of Fig.^ Since this splitting has to overcome 
the Cg bandwidth for a polaron localization, U' is actually 
needed — without U' no strong lattice distortion emerges 
(dotted line), (ii) Due to strong quasiparticle renormal- 
izations the effective bandwidth of the correlated bands is 
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FIG. 2: (Color online) PM spectrum as a function of tempera- 
ture for <; = 0.1 and n = 0.8. With decreasing T, the polaron 
formation is stabilized and a pseudo-gap opens at Ef — 0. 
The dotted line is the PM spectrum without electron phonon 
coupling at n = 0.8, /3 — 30, showing a damped and renor- 
malized quasiparticle peak at Ep- Inset: total FM spectrum 
in the presence (light /green) and absence (dark/red) of lat- 
tice distortion, analogous to the inset of Fig. for the PM 
phase. Due to the enhanced mobility in the FM phase, the 
Jahn- Teller distorted band (i.e., the lower light/green band) 
is wider, not fully occupied, and hence metallic-like. 
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FIG. 3: (Color online) Optical conductivity a-{uj) for the PM 
phase at g = 0.1, p — 30 and n = 0.8, 0.5, and 0.3 electrons/site, 
showing an insulating-like behavior with a reduced weight be- 
low leV. The dotted line shows a metallic Drude peak in the 
absence of the electron-phonon coupling (/3 = 30, n = 0.8); 
the dashed line the bad metallic behavior of the FM phase 
(g = 0.1, /3 = 30, n = 0.8). Inset: The PM resistivity strongly 
increases with decreasing T (here for g = 0.1 and n = 0.8) so 
that a transition to the FM phase reduces the resistivity by 
a factor of ~8. There is a "colossal" magnetoresistance. 



strongly reduced (see the dotted line in Fig. 12). We can 
estimate this bandwidth renormalization to be roughly 
Z = {1 — 9ReE(ijj)/9w|ij=o)~^ ~ 0.4, as obtained from 
the self energy at real lo. Note, however, that these quasi- 
particles are strongly damped through scattering at t2g 
spins (lml](0) ^ — 0.7eV). Since the correlated electrons 
are already more localized, it is easy to trap them as 
polarons. 

Spectrum in the PM phase. Let us turn to the total 
(distortion-averaged) PM spectrum in Fig. [3 It shows, 
upon decreasing T, the development of a pseudogap at 
Ep which separates localized Jahn- Teller polarons and 
undistorted configurations without electrons. This pseu- 
dogap is more enhanced at lower T since — without ther- 
mal smearing of the lattice distortion — the three peaks in 
Fig- n become sharper. Consequently, there is a sharper 
separation (pseudogap) between polaron trapped elec- 
trons and undistorted unoccupied states. 

In Fig. 121 we see that the pseudogap of the PM spec- 
trum leads to a similar pseudogap in the optical spectrum 
at small frequencies, so that the optical spectral weight is 
very much reduced below leV j26i|. The first optical peak 
at ~ 1.5eV (named midgap state in the literature) then 
stems from exciting the polarons into the unoccupied, 
undistorted band above the Fermi energy. At ~ 5— 6eV, 
a second peak in the optical spectrum signals the exci- 
tations of states with two eg electrons per site, which 
correspond to the upper (light/green) peak in the inset 
of Fig. n The reduced optical spectral weight at low 
frequencies, and (approximately) the position of the two 
optical peaks are in agreement with experiment |l2l | , as is 



the shift of the two peaks towards lower frequencies with 
increasing hole doping x (decreasing n= 1 — x) and the 
gradual filling of the optical gap with increasing x [l^ . 
Let us note that the height of the two optical peaks does 
not fully agree with experiment: In contrast to Fig. |31 
the second peak has experimentally several times more 
spectral weight than the first (midgap) peak. We can un- 
derstand this disagreement since there are other optical 
contributions in this energy range: from oxygen charge 
transfer, from La bands, and the upper Hubbard band of 
the t2g electrons. Such effects are beyond our low-energy 
Hamiltonian Also the first (midgap) peak has ex- 
perimentally a larger a{uj) than in Fig. O for a smaller 
(7, would increase and the two peaks would shift to 
lower oj. 

From cr(0), we directly obtain the resistivity in the 
inset of Fig.|31 As the electrons are more strongly trapped 
as Jahn- Teller polarons with decreasing temperature, the 
resistivity of the PM phase is strongly enhanced. We 
observe an insulating-like behavior. 

(Bad) metallic FM phase. Let us finally turn to the 
FM phase which we simulate by fixing the direction of 
the t2g spins. The FM spectrum for n = 0.8 and /? = 30 
is presented in the inset of Fig. |21 separated into con- 
tributions from small and large Jahn- Teller distortion. 
While the peak positions of the distorted and undistorted 
bands are at similar energies as for the PM in Fig. ^ 
the widths of these bands are much larger so that the 
band with Jahn- Teller distortion crosses the Fermi level 
and is hence metallic-like. This is also reflected in the 
self-energy which, in contrast to the PM phase, behaves 
metallic-like with dKeYi{ijj) / doj < 0, corresponding to a 
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quasiparticle weight of roughly However, these 

quasiparticle excitations are very strongly damped since, 
when moving through the crystal, they scatter at differ- 
ent local potentials, with and without Jahn- Teller distor- 
tion. The scattering rate is given by lml](0) ^ — 1.9eV 
which corresponds to a quasiparticle life time of half a 
femto second. It is hence questionable whether we can 
still speak of quasiparticles at all. In any case, a bad 
metallic behavior is to be expected and indeed seen in 
the optical conductivity 'l2] , showing a wide peak at low 
frequencies (dashed line in Fig. ^ . 

Having described the insulating-like PM with a pseu- 
dogap in the spectrum and the FM without this pseu- 
dogap, the magnetoresistance simply follows from the 
PM-FM transition — triggered by temperature or mag- 
netic fields. For n — 0.8, the inset of Fig. O yields a 
change of resistivity by about a factor 8 when going from 
the PM to the FM phase. 

Conclusion. We reported how the tendencies of 
Coulomb repulsion and Jahn- Teller distortion to localize 
electrons mutually enhance each other so that correlated 
electrons are trapped as polarons in the PM phase of 
doped manganites. In the FM phase, on the other hand, 
electrons are still mobile, but strongly scattered because 
of fluctuating (strong and weak) Jahn- Teller distortions, 
yielding a bad metal. The transition from PM to FM is 
hence accompanied by a "colossal" magnetoresistance. 
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